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1 Introduction 

The statistical inference for diffusion models has been thoroughly studied by now; 
see the books [4, 6-8, 10] and references therein. 

In this paper, we consider the homogeneous diffusion process given by the stochas¬ 
tic differential equation 


dX t = 6a{X t ) dt + b(X t ) dW t , 

where W f is a standard Wiener process, and 9 is an unknown parameter. 

The standard maximum likelihood estimator for the parameter 9 constructed by 
the observations of X on the interval [0, T} has the form 
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Ox 


rT a(X t ) 
Jo b(Xt ) 2 
rT g(A't) 

Jo b(X t ) 


dX t 

2 " 

t dt 


see, for instance, [7, Example 1.37] and [9]. If the equation has a weak solution, 

i 2 2 

the coefficient a is not identically zero, and the functions |j-, are locally in¬ 
tegrate, then this estimator is strongly consistent [9, Thm. 3.3]. Moreover, if the 
model is ergodic, then this estimator is asymptotically normal [7, Ex. 1.37]. Note that 
in the nonergodic case the maximum likelihood estimator Or may have different limit 
distributions; some examples can be found in [7, Sect. 3.5]. 

If the data are the observations of the trajectory {X t ,t. > 0} at discrete time 
moments 1 1 , 1 2 ,..., we obtain the discrete-time version of the model. Parameter esti¬ 
mation in such models has been studied since the mid-1980s; see [2, 3, 1 1], A review 
of this problem and many references can be found in [5] and [13]. For recent results, 
see [ 6 , 9, 12], 

In this paper, we are interested in the scheme of observations that is called “rapidly 
increasing experimental design.” The process X is observed at time moments 
ti = iA„, i = 0,..., n, such that A n —»• 0 and nA n —»• 00 as n —» 00 . One of 
possible approaches to parameter estimation is to consider a discretized version of 
the continuous-time MLE Ox- The most general results in this direction were ob¬ 
tained by Yoshida [14]. He proved the consistency and asymptotic normality of the 
discretized MLE in the model, where the process was multidimensional, the drift co¬ 
efficient depended on 0 nonlinearly, and the diffusion coefficient also contained an 
unknown parameter. 

Assume that we observe the process X at discrete time moments = k/n, 
0 < k < n 1+a , where 0 < a < 5 . In this scheme, Mishura [9] proposed the follow¬ 
ing discretized version of the maximum likelihood estimator: 


£fc=o“ a(X±)(Xk±i - X±)/b(X±) 2 
n_1 EfclV a(X±) 2 /b(X±) 2 


She proved its strong consistency in the case where the coefficients a and b are 
bounded. The aim of this paper is to establish the asymptotic normality of this es¬ 
timator. Additionally, we assume the ergodicity of the model, but the boundedness of 
the coefficients is not required. In comparison with general results of Yoshida [14], 
our assumptions are less restrictive. We assume the polynomial growth of the function 
1/6 instead of the condition inf rr b{x) 2 > 0. Also, we do not assume the smoothness 
of the coefficients and the polynomial growth of their derivatives; any Lipschitz con¬ 
tinuous a(x) and b(x) are possible. 

The paper is organized as follows. In Section 2, we describe the model and formu¬ 
late the results. In Section 3, some simulation experiments are considered. The proof 
of the main theorem is given in Section 4. 


2 Model description and main result 

Let (12, jj) be a measurable space. Assume that 0 £ R is fixed but unknown. Consider 
a probability measure Pg such that ]v is Pg-complete. 
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Let X solve the equation 

Xf = xq + 8 f a(X s ) ds + 

Jo 

where xq G R, a, b: M —> M are measurable functions, and { Wt, t > 0} is a standard 
Wiener process on (17,3, Pe). 

Denote c(x) = fj$r,d(x) = jHffa,<Pe(x) = exp{-20 f* c(y) dy},<md<Pg(x) = 
fo <Pe(y) dy. 

Assume that the following conditions hold. 

(Al) For some L > 0 and for any ijgl, 

|a(z) - a(y)\ + \b(x) - b{y) \ < L\x - y\. 

(A2) <Pg(+oo) = -@g(— oo) = +oo. 

(A3) G e := /+“ b{x ^ g ( x ) < °°- 

It is well known that under assumption (Al) the stochastic differential equation 
(1) has a unique strong solution. This assumption also yields that the functions a(x) 
and b(x) satisfy the linear growth condition, that is, 

|a(a;)| + |6(a:)| < Mi(l + |x|) (2) 

for some M > 0 and for all x G R. 

Assume additionally that 

(A4) There exist K > 0 and p > 0 such that 

|&(a?) | 1 < K (l + |a;| p ). 


b(X s )dW s , 


( 1 ) 


Then, for some M 2 > 0 and for any x G K, 

|c(x)| < M 2 (l + |a;| 2p+1 ), \d(x)\ < M 2 (l + \x\ 2p+2 ). (3) 


Under assumptions (A2)-(A3), the diffusion process X is positive recurrent; see, 
for example, [7, Prop. 1.15]. In this case, it has ergodic properties with the invariant 
density given by 


y g (x) = 


X G 


Ggb(x) 2 ipe(x)' 

Let G denote a random variable with density fig(x). Then, for any measurable func¬ 
tion h such that Eg|/i(^)| < oo, 


-y pT /*+oo 

— / h(X t )dt—t / h(x)yg(x) dx = ~Egh(^g) a.s. as T —>• oo, (4) 

' Jo 7-00 


see [7, Thm. 1.16]. Moreover, according to [1, Sect. 11.37], the convergence (4) holds 
also in L\, that is, 

—Eg f h(X t ) dt —> 'Egh(^g). 

1 Jo 

Assume that the invariant distribution satisfies the condition 
(A5) E 0 1^6>r = \x\ r y e {x) dx < oo for all r > 0. 


( 5 ) 
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Let 0 < a < 1. Suppose that we observe the process X at discrete time moments 
= k/n, 0 < k < n 1+a . Consider the estimator 


ZCi c(X^i)AXf 

n-'ECi d{X^)' 


where AX” = X k - XV-i. 

Assume also that 

(A6) a is not identically zero. 

Then E e d(£ e ) > 0. Note also that by (3) and (A5), E gd(^g) < oo. Now we are 
ready to formulate the main result. 

Theorem 2.1. Assume that conditions (A1)-(A6) hold. Then 

(i) 0 n —% 0 as n ^ oo, 

(ii) n a ' 2 (0 n -G)=> N{0 ,1/E fl d(&)) as n —y oo. 

The proof is given in Section 4. 

The following result gives sufficient conditions for consistency and asymptotic 
normality in the case where the parameter 0 is positive. 

Corollary 2.2. Let 0 > 0. Assume that conditions (Al), (A4) and (A6) are fulfilled 
and, additionally, 

limsup c{x) sgn(x) < 0. (6) 

| X | —^OO 

Then 

(i) 0 n -^4 0 as n —> oo, 

(ii) n a ^ 2 {0 n — 0) =>■ N( 0,1/E ed(£e)) asn —>• oo. 

Proof. Note that condition (6), together with (A4), implies that assumptions (A2)- 
(A3) are satisfied and, moreover, all polynomial moments of the invariant density are 
finite; see [7, p. 3]. Hence, the result follows directly from Theorem 2.1. □ 

If the coefficients are bounded, then the consistency and asymptotic normality of 
0 n can be obtained without assumption (A5). 

Corollary 2.3. Assume that conditions (A1)-(A3) are satisfied, the coefficients a(x ) 
and b{x) are bounded, and inf x6 R |6(a;)| > 0. Then 

(i) 0 n 0 as n —> oo, 

(ii) n a ^ 2 (0 n — 0) => N( 0,1/Eed(£e)) as n —>■ oo. 

Sketch of proof. This result can be proved similarly to Theorem 2.1 using the bound¬ 
edness of a(x), b(x), c(x), d(x) instead of the growth conditions (2), (3), and (A4) to- 
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gether with the boundedness of moments of the invariant density. In this case, ( 8 ) im¬ 
plies the inequality 

E e (X t - ) 2m < C(m, 9)n~ m 

n 

for all m £ N and t £ !-, —], k = 1,2,... n a . This estimate is used in the proof 

instead of Lemmas 4.1-4.2. □ 

Remark 2.4. For a £ (0, |), Mishura [9] obtained the a.s. convergence in Corol¬ 
lary 2.3(i) without assumptions (A2)-(A3). 

3 Some simulation results 

In this section, we illustrate quality of the estimator by simulation experiments. We 
consider the diffusion process (1) with drift parameter 9 = 2 and initial value xo = 1 
in three following cases: 

( 1 ) a(x) = 1 — x, b(x) = 2 + sin x, 

(2) a(x) = — arctanx, b(x) = 1 , 

(3) a(x) = -j^r, b(x) = 1. 

Using the Milstein method, we simulate 100 sample paths of each process and find 
the estimate 9 n for different values of n and a. The average values of 9 n and the 
corresponding standard deviations are presented in Tables 1-3. 


Table 1. a(x) = 1 — x, b(x) = 2 + sin a; 




n 






50 

100 

500 

1000 

2000 

5000 

a = 0.1 

Mean 

3.05812 

2.97626 

2.73973 

2.58453 

2.55888 

2.53879 


Std. dev. 

2.06388 

2.00007 

1.43273 

1.34689 

1.26920 

1.22077 

a = 0.5 

Mean 

2.11065 

2.15066 

2.08157 

2.05626 

2.03686 

2.03479 


Std. dev. 

0.62613 

0.56038 

0.31621 

0.28909 

0.22875 

0.18187 

a = 0.9 

Mean 

2.02509 

2.01702 

2.02024 

2.01308 

2.00626 

2.00289 


Std. dev. 

0.27874 

0.19589 

0.09995 

0.06918 

0.04850 

0.03028 




Table 2. 

a(x) = — 

arctan x, b(x) = 1 





n 






50 

100 

500 

1000 

2000 

5000 

a = 0.1 

Mean 

2.69321 

2.66637 

2.65053 

2.66356 

2.59903 

2.46685 


Std. dev. 

2.03142 

2.06075 

1.82903 

1.73034 

1.68212 

1.50186 

a = 0.5 

Mean 

2.12190 

2.10459 

2.01048 

1.99535 

2.01712 

1.99517 


Std. dev. 

0.85304 

0.69484 

0.48803 

0.37807 

0.31746 

0.25846 

a = 0.9 

Mean 

1.95538 

1.97446 

1.98035 

1.99565 

2.00266 

2.00290 


Std. dev. 

0.35057 

0.26796 

0.12235 

0.09050 

0.06496 

0.04533 
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Table 3. a(x) = — » b(x) = 1 




n 






50 

100 

500 

1000 

2000 

5000 

0 = 0.1 

Mean 

1.99507 

1.99813 

1.97122 

1.99255 

1.98366 

1.94811 


Std. dev. 

2.44248 

2.53060 

2.17322 

2.13403 

2.05527 

1.80128 

a = 0.5 

Mean 

1.87038 

1.87897 

1.89022 

1.92593 

1.94964 

1.96624 


Std. dev. 

1.01932 

0.89315 

0.54811 

0.49005 

0.41787 

0.33855 

0 = 0.9 

Mean 

1.90341 

1.92162 

2.00240 

2.00068 

2.00491 

1.99347 


Std. dev. 

0.47656 

0.33693 

0.18136 

0.13173 

0.09595 

0.07033 


4 Proof of Theorem 2.1 


In this section, we prove the main theorem and some auxiliary lemmas. In what fol¬ 
lows, C, C\ , C 2 , • ■ • are positive generic constants that may vary from line to line. If 
they depend on some arguments, we will write C(9 ), C(m, 9), and so on. 

By (1), 

k k_ 

AX%=6 / a{X t )dt + /" b{X t )dW t 


= 9a(Xk~i)— +9 I {aiX t ) - a(Xk- 1 )) dt + b(Xk-i)AWE 

n n J h=± n n 

k_ 

+ f" (b(X t ) — b(Xk~i)) dW t . 


Therefore, 


= 9 + ( c(Xk- 1 )6 / (a(Xt) — a(Xk~ 1 )) dt + 


a(X k- 1 ) 


AWF 


fc=1 


b(Xk=i) ' k 

+ c(Xk- 1 ) / ^ (6(X t ) - &(X*-i)) dWi j / ( ^ ) • 


fc=i 


Then 


9n- 9 = 


“(^ra + + E n ) 

a! : 


(7) 


D n = n _1_ “ ^ 


fc=l 


^ c(Xfc-i )9 r (a(X t ) - a(Xk-i)) dt, 


k=1 
,1+0 


___ a(Xfc-i) 


where 
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E n =J2 c(A^i) / 1 (b(X t ) - b(Xk-i)) dW t . 

h =1 ” ^7T 

Lemma 4.1. Let assumptions (A1)-(A3) and (A5) be fulfilled. Then for every m £ N, 
there exists a constant C(m , 6) > 0 such that 

F,g(X t — Xk-i) 2m < C(m, 9)n~ m+1+a . 

for all n £ N, 1 < k < n 1+a , and t. £ [-AA, ^]. 

Proof. By (1) and the inequality (a + b) 2m < 2 2 m_ 1 (a 2m + b 2m ), 

Ee(A t — Xk~i) 2m 


2 m 


2 m\ 


< 2 


2m— 1 I /)2m 


0 Eg 


z(A s ) ds 


E e | 


b{X s ) dW s 


Using the Burkholder-Davis-Gundy inequality, we obtain 
Eg(X t — A'fc-i ) 2m 


< 2 


2m—1 I r\2m 


Ee 


t \ 2m / t \ m\ 

j ki a(X s ) ds I +C{m)E e [j ki b{X s ) 2 ds\ j. 


By Jensen’s inequality, 


< 2 2m_1 /3^m (+ k—l\ t 


^( 4 _e) 2m X Ee J* i a(A s ) 2m ds 

+ C(m)(t - ^) m_1 Ee J* k b(X s ) 2m ds. 


( 8 ) 


Further, we have 

Ee(A* - Xk~i) 2m < 2 2m_1 ^ 2 m n 1 _ 2 m E 9 a{X s ) 2m ds 

j-n a 

+ C(m)n 1 ~ m Ee / b{X s ) 2m ds 
do 

Now it remains to note that by (2) and (5) the integrals n - “ J” a(X s ) 2m ds and 
n ~ a Jo b(X s ) 2m ds have bounded expectations. □ 

Lemma 4.2. Under assumption (Al), for every in £ N, there exists a constant 
C{m , 0) > 0 such that 

Ee(A t — Afc-i) 2m < C(m,e)n- m E 0 (l + |A^i| 2m ) 
for all n £ N, 1 < k < n 1+a , and t £ — 1 . 

j > — — ’ L n ’ nJ 
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Proof. By ( 8 ), 

E e (X t -X^i) 2m 


<C 1 (m,6)n 1 - m (^ Eg ^ ^a(X s ) 2m ds+ E 9 ^ i ^ s ) 2 m rfsj. 


Using assumption (Al) and (2), we get 


Eg 


/»£ /»£ 

/ i a(X s ) 2m ds < 2 2 m ~ 1 E e / i ((a(X s )-a(X^i)) 2m + a(X*^ 1 ) 2m )ds 

< 2 2 m " 1 LE e (X s - Xfc-i ) 2m ds 

+ 2 2m ~ 1 M(t — -^i)Ee(l + |Xfe-i | 2m ). 

The same estimate holds for Eg fj Li b(X s ) 2m ds. Therefore, 

E fl (X t - X^) 2m < C 2 (m, 9)n 1 ~ m i E e (X s - X*-i) 2m ds 
+ C 2 (m, 0)n- m E e (1 + |X*_i | 2m ), 

and the result follows from the Gronwall lemma. □ 

Lemma 4.3. Assume that conditions (A1)-(A3) and (A5) are fulfilled. Then for any 
m > 1, 1 < * < 2 m, and 0 < j < 2m, dzere exists C(m, 9) > 0 szzc/z that 


Y L 1 E e(|X^i -XtflXtl^dt < C(m,Vw. 

fc=i 

Proof. Applying the Holder inequality and Lemma 4.1, we get 

Ee(|Xfc-i — X t |*|X t | J ) < (E 9 \X^ ~ X t \ 2m+2 )^ (E e \X t \i^^)^^ 

(m-a)i . j(2m+2) . 2m + 2-i 

< C , i(TO,6»)n'T5s+^(E 0 |X t |^+^r) 2m + 2 . 

Then 


Y i_ ± E e(|x^i -x t \'\x t y)dt 


k= 1 


<Cl(m, 0 )n ( 2^+2 / (E^|Xt| 2 m+ 2 ”-i ) 

By Jensen’s inequality we have 


j(2m + 2) . 2m + 2-i 

- 2 ™+ 2 dt. 


. j(2m+2) . 2m + 2-i 

(E e |X t |^+^) 2m+2 dt<n a n““ 


j (2m+2) 

Ee|Xt| 2m + 2-i dt 


2m+2 
2m + 2-z 


By (5) the expression in brackets is bounded. This completes the proof. 


□ 
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Lemma 4.4. Assume that conditions (A1)-(A3) and (A5) are fulfilled. Then 
(i) for any m = 0,1,2 ,..., 


-1 — q \ ' y2m 2/i >2m 

-A-k-i - > E e £ e 


as n —> oo; 


k =1 

(ii) if, additionally, (A4) holds, then 

Li 


D n —4 Egd(^) as n -A- oo. 

Proof, (i) In the case m = 0, the result is trivial. Let m > 1. By (5) we have 

n~ a / X 2m dt ^4 E e $ m as n -> oo. 

Jo 


Hence, it suffices to prove that 


X 2m dt - n _1_a J2 X ^- 


fc=l 


= n~ a E g 


A+" .A 


E lj x l m - x2 &) 


dt 


fe=l 


converges to zero as n -A- oo. By the inequality |x| < \x — y\ + \y\, 


2m 


-y 2 m \<\x- y \ m 2 ™- 1 -* <Y, C >\ x -y\ i \y\ 2 


Therefore, 


t=o 


A+“ 


i =1 

and, by Lemma 4.3, 


F n < Y, °i n ~ a L_ x E e(\ X ^ - XtlW™-*) dt, 


k =1 


F n < C{m.,0)'^2 l Ci n ^+ 2 ’ 0 as n —> 00 . 


i=l 


(ii) For arbitrary x and y. 


d(x) - d(y) = 


a{x ) 2 a(y) 2 


b(x) 2 b(y ) 2 


= (a{x) - a{y)) 

- i b (x) - b{y)) 


i(x) a{y) 


b(x) 2 b(x)b{y) 
a{x)a{y) a{x)a(y) 
b(x) 2 b{y) b(x)b(y) 2 J 


(9) 


By (Al), (A4), and (2), 
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| d(x) - d(y )| < C|a; - y |(l + |a ;| 2p+1 + (l + M P )(l + \y\ p+1 ) 

+ (i + M 2p+1 ) (i + \y\ p+1 ) + (i + M p+1 ) (i + |y| 2p+1 )) ■ 

(10) 

The rest of the proof can be done similarly to part (i) using estimate (10) instead 
of (9). □ 

Lemma 4.5. Under the assumptions of Theorem 2.1, 

(i) n~ a / 2 \A n \ 0 as n —> oo, 

(ii) nr a / 2 \E n \ 0 as n —> oo. 

Proof, (i) By the Cauchy-Schwarz inequality we have 

n 1+ “ „A 

E 0 K|<|0|^ fE,! c{Xk-i)(a{X t ) - a(Xfc-i)) 

fc=l 

■n}+ a .A 


dt 


<|0|£) / i (E@c(Xfc-i ) 2 ) 5 (Eg (a(Xt) — a(Xfc-i)) 2 ) 5 dt. 


k= 1 

Using (Al), (3), and Lemma 4.2, we get 

n 1+ “ .A 


E e |A n |<C 1 |0| ^ / i (E e (l + |X^- i | 4 J ,+2 )) 5 (E 0 (Xt-X^- i ) 2 ) 5 df 


fc=l 




< C 2 (9)n~ 1 / 2 Y / (E fl (l + |X*__i| 4p+4 )) 


1/2 


dt 


fc=l 
,l+o 


= C 2 (0)rr 3 / 2 £ (E e (l + |*^i| 4p+4 )) 


1/2 


fc=l 


< C 2 (9)n a ~ 1 ^ 2 ( n _1_a ^ E e (l + |AA-i | 4p+4 ) 


1/2 


fc=i 


By Lemma 4.4 the expression n 1 a Ylk=i Eg(l + |A'fc-i | 4p+4 ) is bounded. There¬ 
fore, 

rr“ /2 E 0 |A„| < C 3 (0)n5( Q:_1 l ^ 0 as n -+ oo. 


(ii) We have 


E n = f h(t, ut) dWt, 
Jo 


where 
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Then 


h(t,Lu) 2 dt = /“ E 0 (c(X^) 2 (6(X t )-6(X^i)) 2 )dT 

"'° fe=i 

Similarly to (i), we can estimate E gE 2 < C^{d)n a ~ 1 . Therefore, n~ a J 2 \E n \ 0 
as n —> oo. □ 

Lemma 4.6. Under the assumptions of Theorem 2.1, 

p 


(i) n a B n —0 as n —> oo, 

(ii) n~ a l 2 B n => N( 0, E gd(£g)) as n —»• oo. 

Proof, (i) Let us prove the convergence in L 2 . We have 


where 
B' = 


a(X t ) 


Jo b{X t ) 

Then by (5) we have 


B n — B' n + B ", 


dwb " =e /:, 


a(X t ) 


-[Jh^KKX^) b(X t ) 


dW t . 


C r 

E e (n- a / 2 B' n ) 2 =n~ a Eg / 

Jo 


a{X t y 
b(X t ) 2 


f-n a 

dt = n~ a E 6 / d(X t ) dt’Egditg) 
Jo 

( 11 ) 


as n —> oo. Hence, n a B' n 0 as n —>• oo. 

Arguing as in the proof of Lemma 4.5 (ii), we obtain 


E «K ) 2 = £ / , 


fc=l 


a(X t ) 


b(X^i) b(X t ) 


dt. 


Further, 


a(x) a(y) a{x) , , 1 , , , , 

W)~m = wm m _ 6 W)+ m' (l (} ~ (s)) ' 

Therefore, by (2) and assumption (A4), 

a(x) a(y) „ 2a(:r) 2 . ,, ,>2 2 , ,2 

< C(a: - y) 2 (l + |x| 2p+2 |t/| 2p + |t/| 2p ). 

Similarly to the proof of Lemma 4.4, we get the convergence 

n _ “Eg(H ") 2 —» 0 as n —>• 00 . 

(ii) According to [7, Theorem 1.19], it follows from (11) that n~ a ^ 2 B' n => iV(0, 
Eed(^e)) as n —> 00 . Taking into account the convergence (12), we obtain the result. 

□ 


( 12 ) 


Now the statement of Theorem 2.1 follows from (7) and Lemmas 4.4-4.6. 















28 


K. Ralchenko 


Acknowledgments 

The author is grateful to the anonymous referee for his careful reading of this paper 

and suggesting a number of improvements. 

References 

[1] Borodin, A.N., Salminen, P.: Handbook of Brownian Motion: Facts and Formulae, 2nd 
edn. Birkhauser, Basel (2002) MR1912205. doi:10.1007/978-3-0348-8163-0 

[2] Dacunha-Castelle, D., Florens-Zmirou, D.: Estimation of the coefficients of a diffusion 
from discrete observations. Stochastics 19, 263-284 (1986) MR0872464. doi:10. 1080/ 
17442508608833428 

[3] Florens-Zmirou, D.: Approximate discrete-time schemes for statistics of diffusion pro¬ 
cesses. Statistics 20(4), 547-557 (1989) MR1047222. doi:10.1080/02331888908802205 

[4] Heyde, C.C.: Quasi-Likelihood and Its Application. A General Approach to Optimal Pa¬ 
rameter Estimation. Springer, New York, NY (1997) MR1461808. doi:10.1007/b98823 

[5] Iacus, S.M.: Simulation and Inference for Stochastic Differential Equations. With R Ex¬ 
amples. Springer, New York, NY (2008) MR2410254. doi: 10.1007/978-0-387-75839-8 

[6] Kessler, M., Lindner, A., Sprensen, M. (eds.): Statistical Methods for Stochastic Differen¬ 
tial Equations. Selected Papers Based on the Presentations at the 7th Seminaire Europeen 
de Statistiques on “Statistics for Stochastic Differential Equations Models,” La Manga 
del Mar Menor, Cartagena, Spain, 7-12 May 2007. CRC Press, Boca Raton, FL (2012) 
MR2975799 

[7] Kutoyants, Y.A.: Statistical Inference for Ergodic Diffusion Processes. Springer, London 
(2004) MR2144185. doi:10.1007/978-l-4471-3866-2 

[8] Liptser, R.S., Shiryayev, A.N.: Statistics of Random Processes. II. Applications. 
Springer-Verlag, New York, Heidelberg, Berlin (1978) MR0488267 

[9] Mishura, Y.: Standard maximum likelihood drift parameter estimator in the homoge¬ 
neous diffusion model is always strongly consistent. Stat. Probab. Lett. 86 , 24—29 (2014) 
MR3162713. doi: 10.1016/j.spl.2013.12.004 

[10] Prakasa Rao, B.L.S.: Asymptotic Theory of Statistical Inference. lohn Wiley & Sons, 
New York (1987) MR0874342 

[11] Prakasa Rao, B.L.S.: Statistical inference from sampled data for stochastic processes. 
Contemp. Math. 20, 249-284 (1988) MR0999016. doi:10.1090/conm/080/999016 

[12] Sprensen, M.: Parametric inference for discretely sampled stochastic differential equa¬ 
tions. In: Handbook of Financial Time Series. With a Foreword by Robert Engle, pp. 531— 
553. Springer, Berlin (2009) 

[13] Sprensen, H.: Estimation of diffusion parameters for discretely observed diffusion pro¬ 
cesses. Bernoulli 8(4), 491-508 (2002) MR1914700 

[14] Yoshida, N.: Estimation for diffusion processes from discrete observation. I. Multivar. 
Anal. 41(2), 220-242 (1992) MR1172898. doi:10.1016/0047-259X(92)90068-Q 


